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Abstract 

The dual parametrization of generalized parton distributions (GPDs) is con- 
sidered in details. We discuss which part of information about hadron structure 
encoded in GPDs [part of total GPD image] can be restored from the known am- 
plitude of a hard exclusive process. The physics content of this partial image is 
analyzed. 



Introduction and basics of the dual parametrization of GPDs 

The generalized parton distributions (GPDs) [1] (for recent reviews of GPDs see 
Refs. [2, 3, 4]) describe the response of the target hadron to the well-defined QCD quark 
and gluon probes of the type: 
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where the operators are defined on the light-cone, i.e. z 2 = 0, and A, B are various 
hadronic states. In this way the hard exclusive processes provide us with a set of new 
fundamental probes of the hadronic structure. The hadron image seen by the probes (1) 
is encoded in the dependence of GPDs on its variables. 

The amplitudes of hard exclusive processes, which we consider as direct observables, 
are given by the convolution of GPDs with perturbative kernels. It means that mea- 
surements of the amplitudes provide us with some kind of sectional images of GPDs - 
the convolution integral "projects out" one of variables in GPDs. In this contribution 
we address the question: What part of the full hadron image provided by GPDs can be 
reconstructed if we know the amplitude of hard exclusive process? This question was 
recently discussed in Refs. [5, 6, 7, 8]. Here we follow the Ref. [7] adding new ideas and 
calculations. 

For our analysis we employ dual parametrization of GPDs suggested in Ref. [9]. This 
parametrization is based on representation of parton distributions as an infinite series of 
t-channel exchanges [10]. 
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In the dual parametrization the GPD H(x, £, t) is expressed in terms of set of functions 
Qivip, t) [y = 0, 1, 2, . . .). For the detailed expression see original paper [9]. We call the 
functions Q 2u (x, t) forward-like because: 

• At the LO scale dependence of functions Q 2l/ (x, t) is given by the standard DGLAP 
evolution equation, so that these functions behave as usual parton distributions 
under QCD evolution. 

• The function Q (x,t) at t = is related to the forward distribution q(x) as: 



X f ClZ 

Q (x,t = 0)= (q(x) + q(x)) - - — (q(z) + q{z)) 

L Jx Z 

The inverse transformation reads: 

q(x) + q(x) = Q (x) + ^ £ A Q o(y ) . 



(2) 



(3) 



• The expansion of the GPD H(x,£,t) around the point £ = with fixed x to the 
order £ 2u involves only finite number of functions Q 2fM (x, t) with /i < v. For instance, 
to the order £ 2 these are only Q and Q 2 : 



H(x, Qo(x, t) + ^ £ A g o(y , t) 
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The Eq. (4) reveals some of important properties of the forward-like functions. For exam- 
ple, Eq. (4) dictates that the forward like functions Q 2u (x, t) can have small x singularities 
of the type Q 2u (x,t) 

^ 2v+a where the exponent ot governs small x behaviour of forward 
quark distribution q(x) + q(x) ~ l/x a . At first glance such singularities lead to divergency 
of of the D-term and of the Mellin moments of GPDs. However one can easily show that 
the most singular at small x terms of Q 2u (x,t) has the following form: 
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f 1 dz z k R u {z) =0, k<2v-l. 
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Here we introduce factor containing regularizing parameter e which is put to zero at 
the end. Introducing the regularizing factor we modify the parton distribution in the 
arbitrarily small vicinity of x = 0, after performing calculations one has to check that 
the final result has a finite limit at e — > 0. The specific form (5) of the small-x singular 
terms of forward like functions Q2v(x,t) (with v > 1) is actually the consequence of the 
polynomiality of the GPDs Mellin moments. The above prescription is very important to 
understand the naively divergent integrals in our expressions below. The functions R2v(z) 
can be computed if one knows the expansion of the GPD at small £: 

H(x,£,t)=q(x,t)+e f 2 (x,t) + ... (7) 

For example for the simple model of GPD H(x,£) = q(x) we obtain that: 

R x (z) = z 6'{1 - z) 5{l- z) - [1 + 3z) . (8) 

We shall see below that the amplitudes of the hard exclusive processes are expressed in 
terms of single function 

oo 

N(x,t) = Y,a? v Q 2v (x,t). (9) 
Its small x behaviour has the following form 



tfOM)~^ + ..., (io) 

with a < 2. Important to note that neither the exponent a nor the coefficient No in front of 
l/x a is determined by the forward parton distributions. However, as we discussed above, 
the singular terms, which deviate from those dictated by the forward parton distributions 
t have very specific form (16). In the first phenomenological applications of the dual 
parametrization to describe the DVCS data [11] this was not taken into account (it was 
assumed that Q2v{x) ~ l/x a ). The inclusion of the singular terms of the type (5) into 
the parametrization can considerably change the results and conclusions of these papers. 



Forward-like distributions in terms of amplitudes 

The leading order amplitude of hard exclusive reactions is expressed in terms of the 
following elementary amplitude''': 



dx H(x, £, t) 



Z-x-iO £ + x-i0 



+If the forward distribution behaves q(x) <~ l/x a the corresponding peace has the following behaviour 
N{x,t)~%£ l/x° 

■i-We restrict ourselves to the singlet (even signature) amplitudes, generalization for odd signature 
amplitudes is trivial. 
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We see that the amplitude is given by the convolution integral in which dependence of 
GPDs on variable x is "integrated out". Mathematically from the equations (11) one 
can not completely restore the GPD H(x,£,,t). So we are not able to perform "complete 
imaging" of the target hadron from the knowledge of the amplitude. The key question 
is: what part of the "complete image" can be restored from the known amplitude? What 
is the physics content of the restorable part of the complete image? Attempt to answer 
these questions in the framework of the dual parametrization was done in Ref. [7]. We 
refer the reader to this paper for details, here we just give main results and more detailed 
discussion of subtleties. 

We can express the amplitudes in terms of forward-like function (see Eq. (9)) N(x,t) 
as following [9]: 
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Here wc introduced the D-form factor: 



D(t) = ±d n (t)= 1 - f 1 dz^^, 

71=1 Z J ~ 1 1 Z 



(13) 



where D(z,t) is the D-term [12]§. 

Now we clearly see that the knowledge of the LO amplitude is equivalent to the 
knowledge of the function N(x,t) and D-form factor D(t). Moreover the D-form factor 
can be computed in terms of N(x, t) and Qo(x, t). Note that the latter function is to great 
extend is fixed by the forward parton distributions, see Eq. (2). The expression for the 
D-form factor is the following [7] : 



D{t)= f 1 -Q (z,t) f—L=-l) + f ^ [N(z,t)-Q (z,t)} 
Jo z \yl + z l J Jo z 
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(15) 



Given the small x behaviour (10) of the function N(x,t) naively one may think that the 
above integral is divergent. However it is not the case due to the very specific form of the 
singular terms, see Eq. (5). Because of this specific form of the small x singularities we 
can represent the function N(x, t) in the following way: 



§ The coefficients d n {t) individually can be obtained from the generating function [7]: 

oo 1 f 1 rl °° / 1 \ 



(14) 
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N(x, t) = Q {x, t) + N^(x, t) + N shl %(x, t) , 

where the function N Teg (x,t) goes to zero as x — > whereas the singular part behaves as 
~ l/x a and has the following structure: 

N**(x,t) = £ | Q„(M) ^ i2 g ) , (16) 
where the function R(z) has the property: 

r 1 1 

/ dz- R(z) = 0. (17) 
Jo z 

Now we can rewrite the representation for the D-form factor (15) as follows: 

= £ f [Q„(M) + iV^(M)] (^-1)+/; f **■(*,.) -^.(18) 

Here all integrals are explicitly convergent and therefore one can safely put e = 0. It 
implies that, for practical calculations, one can use the following prescription. If one 
knows the small x behaviour of the function N(x,t) — Q (x,t)^: 

N(x,t) -Q (x) ~E4J + ---> with0<a,<2 (19) 



then we can define the singular part of the function N(x, t) in the following way: 

with arbitrary function R(z) which satisfies Eq. (17). Now it is obvious that the function 
defined as: 

iV reg (x, t) = N(x, t) - Q (x, t) - N si ^(x, t) , (21) 

tends to zero if x — > 0. 

As it is shown in Ref. [7] the Mellin moments of N(x, t) allow to fix the angular 
momentum of exchanged partons. 

^ dx x J ~ l N(x, t) = X - J 1 dz ®^*'^ , (22) 



^We see below that it is equivalent to the knowledge of the small £ behaviour of the amplitude which 
can be directly measured and small x asymptotic of the forward parton distribution 
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where t) is the distribution amplitude corresponding to two quark exchange in the t- 

channel with fixed angular momentum J. The quantity on RHS of Eq. (22) carries valuable 
information about the hadron structure - it tells how the target nucleon responses to the 
well defined string-like quark- ant iquark probe (see Eq. (1)) with fixed angular momentum 
J. 

We spent some time to discuss how to compute the D-form factor if one knows the 
functions N(x, t) and Qo(x, t). The latter function is closely related to the forward parton 
distributions and the former can be computed directly from the observed amplitudes. The 
final result for the inversion problem is the following (see details [7]): 



7T 



This remarkable formula allows to restore the function N(x,t) from the measured imagi- 
nary part of the amplitude. Note that the inversion formula contains the amplitude only 
in the physical region. The Eq. (23) directly translates any knowledge (theoretical, model 
or experimental) of the amplitude into comprehension of GPDs. At £ — > 1 the imaginary 
part of the amplitude should go to zero. Let us assume that Im A(£,t) ~ (1 — as 
£ — > 1, than from Eq. (23) one can easily obtain that 

as x — > 1. 

As to small x behaviour, the Regge like asymptotic of the imaginary part of the 
amplitude Im A(£,t) ~ l/£ a , according to Eq. (23), corresponds to the following small x 
behaviour of N(x,t): 



N(x,t) 



2 a r(|)r(i + «) x a ' 

Up to now we made use of only the imaginary part of the amplitude, let us study the 
real part of the amplitude. For this we take the expression (23) for the N(x, t) in terms 
of ImA(£,i) and substitute it into Eq. (12). After some simple calculations we arrive to 
the following expression for real part of the amplitude: 



ReA(C,t) = 2D(t) + -vp f d( lmA((,t) 

IT JO 

in which we can immediately recognize the dispersion relation for the amplitude with 
one subtraction at non-physical point £ = oo (corresponding to v — (s — u)/Am = 0). 
The D-form factor is the corresponding subtraction constant. This result was obtained 
recently in Refs. [14, 15] by independent methods. We see that the dual parametrization 
automatically ensures the dispersion relations for the amplitudes. 

The very idea of the dual parametrization of GPDs in terms of t-channel exchanges 
was motivated by the crossing relations [10] between GPDs and generalized distribution 
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amplitude [16]. The later enter the description of the hard exclusive processes in the 
cross channel, like 7* + 7 — + h + h. In the LO the amplitude of the cross process can be 
expressed in terms of the function N(x,t) analytically continued to time-like t (t > 0): 

A cross { V ,t) = / —N(x,t) 
Jo X 

Here — rj is directly related to the cos(0 cm )- cosine of scattering angle in centre of mass 
system, see for details Refs. [16, 10]. Now substituting our inversion formula (23) into 
this expression we obtain, rather simple result: 

2 [W £ ( £ \ 

A c ^( v , t) = - dt-±- lmA(±-t)+2 D(t) . (26) 

TV Jo 1-4 VM / 

Actually this equation is the consequence of the dispersion relation (24). 

Now we have all ingredients to make proposals for possible ways to analyze date on 
hard exclusive reactions. 

Possible ways to comprehend GPDs from data 

The observables of hard exclusive processes can be analyzed assuming certain func- 
tional form of lmA(£,t) with set of free parameters. It can be the form with motivated 
by some theories (like Regge-motivated [17], etc.) or one can choose some flexible enough 
functions. For instance one can take as the functional form: 

lmA(Z, t) = Wo(C, *) + £ 4 (! - > ( 27 ) 

i x 

where Pi, fa is the set of free parameters which should be adjusted to data, ImA (£, t) 
is the amplitude computed with the first Eq. (12) taking N(x) = Qq(x), i.e. it is fixed 
by the forward parton distribution. The reader can select any other function form of 
the amplitude. With help of our master formula (23) we can immediately compute the 
function N(x,t), for instance for the choice (27): 

N(x,t)=Q (x,t) + Y i P i F^(x,t), (28) 

i 

with F^ A \x,t) given by the simple combination of hypergeometric functions. Knowing 
the function N(x, t) and the forward parton distribution Qo(x, t) we can calculate the real 
part of the amplitude and the D form factor. To compute the D form factor we should 
figure out the small x asymptotic (19) and construct N sms (x,t) according to Eq. (20). 
For the amplitude of the form (27) we obtain: 



a/1 - 2xi] + x 2 + a/1 + 2x7] + x 2 a/IT 



x- 



+ 2D(t) . (25) 
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N^(xt) = VP + ^ [ 1 Si dz z^R(z) 

y,J Y l 2^V^ r K + i) 1^ fodzz^Riz) 
A(2^-l) 1 f*dz z^Rjz) 1 
ai x 01 *- 1 Jo 1 dzz^- 2 R{z) J ' 



(29) 



with arbitrary function i?(z) possessing the property (17). After one can compute the 
D form factor (18) and the real part of the amplitude (12) in terms of the adjustable 
parameters which enter the chosen functional form for the imaginary part of the amplitude. 
In principle, these adjustable parameters can be determined from the data on imaginary 
part of the amplitude (like beam spin asymmetry for DVCS), then the real part and the 
D- term is predicted. Again let us give the result for our choice of the fitting function 
(27). For such choice the D form factor is: 



m = /; f Qo( ,, t) ( i _ x) + i E i ^^+f oo) 

JO Z Vvl + -2 / 7T i L(l+/Ji-ai) 

Till now we assumed that the forward limit of the GPD at hand is known. However in 
some cases the forward limit is not known and should be determined. The most interesting 
example is the nucleon GPD Hm(x, £, t) = H(x, £, t)+E(x, £, t)K Its forward limit is given 
by unknown function Qq^ (x, t) which is normalized to the angular momentum carried 
by quark in the nucleon J Q : 



j 1 dxxQ i i 1 \x,t) = l J Q (t). (31) 
Jo 3 

Actually it is one of the main goals of measurements of hard exclusive processes to deter- 
mine this function. 

Well, as we do not know the forward function Qq M ^ (x, t) there is no point to single out 
its contribution to the imaginary part of the amplitude. Instead, we simply parametrize 
the Im ^4( M )** as the whole, e.g. by the following function: 



Iim4W(£,t)=£ (32) 

i x 

Now again we can easily perform analysis along the line discussed above, the only differ- 
ence is that the expression for the D-form factor is different: 



D(t) = -f 1 - [Q ( M \z,t)-Qi M >^(z,t)]+^ Pl 

JO Z 1 J 7T Y 



T(-a t ) r(i + A) 
r(i + A-0 



(33) 



"This "magnetic- like" GPD has correct angular momentum quantum numbers in the i-channel. The 
full classification of such nucleon GPD is given in [3] . 

**We keep the superscript M to stress that we consider the case with unknown forward limit 
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Here Qq M \x, t) is unknown (to be determined) forward parton distribution and Q M ^ sms (x, t) 
its singular part computed according to the Eq. (20) with obvious replacement of func- 
tions. Now we can parametrize the unknown function Q M ^ (x, t) and vary its param- 
eters to achieve the agreement with data. In the case of nucleon GPD i/^x, £,t) = 
H(x, £, t) + E(x, £, t) we know that corresponding D-form factor is zero and we can con- 
strain the forward function Q M ^ (x, t) from the following equation: 



I 



J l{l + IJi-ai) 



o z 



Here the parameters Pi, ai and /3j are determined from fitting the imaginary part of the 
amplitude 1 " 1 ". Note that the expression in Eq. (33) corresponds to the D-form factor 
computed from the amplitude (32) with help of analytical regularization discussed in 
Refs. [5, 6, 8]. 

Up to now we discussed possibilities to fit data at fixed photon virtuality Q 2 . If we 
are interested in studies of the scaling violation and determination of the whole GPD we 
have to model the forward-like functions Q2u{x,t) individually. In this case one has more 
freedom and one needs more a priori input. 

For modelling of forward-like function we propose the following general ansatz: 

Q*,(x, t) = f - Q (z, t) P v (-) + 4 f ~Qo (-, t) z 2 »R v (z) . (36) 

J x Z \ Z / 0C J x Z \ Z / 

Here the second terms accounts for most singular small x asymptotic of the forward-like 
functions. As we discussed above, the specific structure of these singularities (5) ensures 
that R v (z) satisfies the condition: 

I -I^k z2UR Az) = ^ k = l,3,...,2i/-l. (37) 

Possible ways to model the functions P v {z) in Eq. (36) were discussed in Ref. [7]. As to 
modelling functions R v {z) one may try the following form: 

2v-\ 

R v {z) =w(z)J2 r j Pj( z ), (38) 

j=0 



^The reader can easily derive corresponding equation for her/his beloved functional form of the am- 
plitude. As a curiosity :) we note that if we assume that forward limit of "magnetic" nucleon GPD has 
the form: 

Qi M \x,t) = ^ a {l- x )\ 

we can obtain J® from the fit to data on lmA^ M \^,t) by the functional form (32) as follows: 

jQ = 3a( g -l) v Tj-Oi) r(l + ft) 

57r(2 + 6-a)(l + 6-a) ^ 1 Y{1 + fa - a,) ' [ ' 
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where Pj(z) is set of orthogonal on the interval [0, 1] with the weight w(z) polynomials. 
This form automatically satisfies the condition (37). Alternatively one can choose to 
model the functions R u (z) allowing the presence of (^-functions, like in Eq. (8). 

We see that there there is rather large freedom in choice of the functional form for 
the forward-like functions Q2v(x). Further constraints for possible form and size of these 
functions we can obtain computing the forward-like functions in various models for GPDs, 
like, for instance, chiral quark-soliton model [19]. 

Further possible uses of dual parametrization 

Here, instead of conclusion, we outline couple of ideas about use of dual parametriza- 
tion apart from those which obviously follow from the discussion above. 

One of advantages of the dual parametrization is that it allows to separate the con- 
tribution of the usual forward parton distributions to the amplitude of a hard exclusive 
process. This is the function Qq(x, t) in the function N(x, t) = Qo(x, t) +x 2 Q2(x, t) + . . ., 
the latter completely determines the amplitude. This feature allows us to study "truly 
non-forward" effects in the hard processes. For example, studying DVCS on nuclei, with 
help of dual parametrization, we can easily single out "non-forward EMC effect". The 
usual EMC effect is taken into account by the function Qo{x, t) which is determined by nu- 
clear forward parton distributions. Model calculations [13, 18] predict unusual behaviour 
of x- moment of the function Q2(x, t) with atomic number. The dual parametrization 
allows to reveal such new nuclear effects from the data. 

Another possible application of the dual parametrization is related to Eq. (22). Apart 
from obvious uses of the Eq. (22), it can bring new insight into spectroscopy of baryon 
resonances. One can repeat the discussion above for the "non-diagonal" GPDs which 
are given by the matrix elements of the light-cone operators (1) for transitions between 
nucleon and meson- nucleon states**, in this case one can introduce the function N(x, t, . . .) 
which in this case depends on additional variables (denoted by . . ., the definition of these 
variables and their relation to partial wave decomposition of the final meson-nucleon 
system can be found in Ref. [20]) which characterize the final meson-nucleon state (its 
invariant mass, orbital momentum, etc.). The corresponding function can be obtained 
from data using the Eq. (23). Now if one computes the x J_1 moment of this function, 
one obtains, according to Eq. (22), the meson-nucleon state produced by the well defined 
probe with spin J. This allows to excite new baryon resonances which couple weakly to 
standard probes like photons, also it allows to study the properties of known resonances 
comparing their excitations by probes with different spins. 
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